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Upper bounds on the bulk energy dissipation rate e and enstrophy dissipation rate \ are derived 
for the statistical steady state of body forced two dimensional turbulence in a periodic domain. For 
a broad class of externally imposed body forces it is shown that e < kfU 3 Re~ x l 2 (C\ + C^-Re^ 1 ) 
and x < ^/^ 3 (Ci + C^-Re -1 ) where U is the root-mean-square velocity, kf is a wavenumber (inverse 
length scale) related with the forcing function, and Re = U /vkf. The positive coefficients Ci and C2 
are uniform in the the kinematic viscosity v, the amplitude of the driving force, and the system size. 
We compare these results with previously obtained bounds for body forces involving only a single 
length scale, or for velocity dependent a constant-energy-flux forces acting at finite wavenumbers. 
Implications of our results are discussed. 

PACS numbers: 



I. INTRODUCTION 

The study of two dimensional turbulence was originally 
justified as a simplified version of the more complex three 
dimensional turbulence, but it has come to be regarded 
as an interesting research field in its own right with deep 
connections to geophysical and astrophysical problems 
such as as strongly rotating stratified flows Q. A large 
number of experimental methods have been devised to 
constrain flows in two dimensions (e.g. soap films) al- 
lowing some theories theories to be tested in the lab 0- 
Direct numerical simulations are far easier than the three 
dimensional case, and this has enabled researchers to in- 
vestigate two dimensional turbulence computationally at 
much higher Reynolds numbers |E El IE IE 0> As a 
result, there are more simulation data to test theories of 
two dimensional turbulence. Nevertheless many funda- 
mental questions remain open; see [| for a recent review. 

The inviscid conservation of enstrophy as well as en- 
ergy in two dimensions results in two cascading quadratic 
invariants that make the phenomenology of two dimen- 
sional turbulence somewhat more complex than three di- 
mensional turbulence and not derivable from simple di- 
mensional arguments. Theoretical studies of turbulence 
usually employ a statistical description and often involve 
assumptions about homogeneity isotropy and the nature 
of the interactions. Based on physical arguments, Kraich- 
nan 0, Leith ^(j and Batchelor E3 conjectured that 
there is a dual cascade in two dimensional turbulence: 
energy flows to the larger scales while enstrophy moves 
to the small scales (when the system is driven at some 
intermediate scale). Kraichnan-Leith-Batchclor (KLB) 
theory assumes isotropy and homogeneity in the limit of 
infinite domain size, and in the zero viscosity limit pre- 
dicts a A: -5 / 3 energy spectrum for the large scales and a 
k~ 3 energy spectrum for the small scales. The assump- 



tions of the KLB theory, as well as the power index and 
the universality of the two energy spectra has been ques- 
tioned in the literature [H [Hill HE EE El El 

In this paper we derive some simple rigorous bounds 
for the long time averaged bulk energy and enstrophy dis- 
sipation rates for two dimensional statistically stationary 
flows sustained by a variety of driving forces. The study 
of physically relevant rigorous bounds on the energy dis- 
sipation rate, i.e., the power consumption of turbulence, 
for a class of boundary-driven flows can be traced back to 
the seminal work of Howard |l9l | , and in recent years an 
alternative approach |20l| has renewed interest in those 
kinds of problems, providing direct connections to exper- 
iments in some cases. Bounds for the energy dissipation 
of steady body forced flows — more convenient for theo- 
retical and numerical investigations — in three dimensions 
have been derived by Foias ^ and others pll23Ll2ll25l|. 
Not unexpectedly, Foias et al. [2(| also derived a bound 
for the enstrophy dissipation rate in the statistically sta- 
tionary states of two-dimensional turbulence driven by 
a restricted class of forces. Bounds for the energy and 
enstrophy dissipation in two dimensional flow driven by 
a monochromatic forcing were derived in |27| and |l6j| . 
The case of temporally white noise forcing was studied 
by Eyink [2^ . More recently Tran el at [22, H(j derived 
bounds of the enstrophy dissipation for freely decaying 
two dimensional turbulence in terms of the initial ideal 
invariants of the flow. (See also [3lJ for the treatment 
of the same problem in terms of the inviscid Euler Equa- 
tions). Finally, we mention that bounds on the dimension 
of attractor of the 2D Navier-Stokes have been derived 
in miH Hi HE HI and more recently by ^ . 

The results for the energy and enstrophy dissipation of 
forced flows derived in this paper apply to a more gen- 
eral type of forcin g th an the single scale forcing that was 
considered in fltl l27l|. We also consider forces that are 
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smoothly varying in time, unlike the white noise forcing 
investigated in J23 , and we are particularly interested in 
the behavior of the long times averaged dissipation rates 
in the vanishing viscosity limit. Because the viscosity is 
a dimensional quantity we must specify what we mean 
by "small" viscosity. To be precise, we measure the mag- 
nitude of the viscosity in terms of the Reynolds number 
in the statistically steady state, 



Re 



U 

k f v 



(1) 



where U is the root mean square velocity and kf is a 
natural wavenumber (inverse length scale) in the driving 
force. The dissipation rates are also dimensional quanti- 
ties, so we measure them in terms of the inviscid scales 
determined by U and kf. That is, we estimate 



ft 



k f U 3 



and 



7 



X 



k 3 f U 3 



(2) 



in terms of Re and focus on the Re — > oo limit holding 
other parameters (such as the large scale eddy turnover 
time in the most general case) fixed. For a broad class of 
external driving we find that 



ft < Re- 1/2 and 7 < Re , 



(3) 



consistent with an enstrophy cascade of sorts. 

However, for special cases of forcing such as "ultra 
narrow band" monochromatic (i.e., involving on a single 
length scale, albeit with a broad range of time depen- 
dence) forcing, or for a fixed energy flux forcing popular 
for direct numerical simulations, a stronger bound holds: 



ft Z Re- 



sold 



7 ^ -Re" 



(4) 



This kind of Re^ 1 scaling suggests "laminar" flows where 
the energy is concentrated at or above the smallest length 
scale of the forcing. This sort of scaling has been shown 
before in the literature [l6ll2^ | for monochromatic forcing 
and for white noise in time forcing psj . 

In every case the bounds derived here are strictly less 
than those available -or expected- for three dimensional 
turbulence. The upper bounds J3J on the energy and 
enstrophy dissipation for two dimensional flows derived 
here are in a sense a consequence of combining the ap- 
proaches in [i^l and 0, |2^| applied to a class of forcing 
functions concentrated in a finite range of length scales. 
Even though some steps in our analysis have been taken 
before, in order to make this paper self-contained the 
complete (but nevertheless short) proofs will be presented 
here. 

The rest of this paper is organized as follows. In the 
next section II we introduce the problem and basic defini- 
tions, and perform the analysis leading to © for the sim- 
plest case of time independent body forces. Section III 
generalizes the analysis to include a broad class of time- 
dependent forces. In section IV we briefly review the 



results for time-dependent but single-length scale forcing 
and "fixed-flux" forces in order to establish the stronger 
results in 10}. The concluding section V contains a brief 
discussion of the results and their implications. 



II. TIME-INDEPENDENT FORCING 

Consider a two dimensional periodic domain [0, L} 2 , 
i.e., T 2 , filled with an incompressible fluid of unit density 
evolving according to the Navier-Stokcs equation: 



<9 f u + u ■ Vu 



-Vp- 



jA7 2 u 



(5) 



where u = iu x (x,y,t) + ju y (x,y,t) is the incompressible 
(divergence-free) velocity field, p(x, y, t) is the pressure, 
v is the viscosity, and f(x, y) = if x (x,y) +jf y (x,y) is 
a smooth, mean zero, divergence-free body force with 
characteristic length scale ~ kf (defined precisely below). 
The scalar vorticity u> = d x u y — d y u x satisfies 



d t uj + u • Vw = vV 2 uj + <fr 

where (j> = k-V x f = d x f y - d y f x . 

The Reynolds number is defined in Q where 

E/^dul 2 ) 1 / 2 



(6) 



(7) 



is the root-mean-square velocity with (•) representing the 
space-time average 

<g> =^ oo ^J q (j? J r2 g(*> y, t) d 2 x^j dt. (8) 

(The limit in the time average is assumed to exist for 
all the quantities of interest.) The forcing length scale 
associated with the wavenumber kf is defined by 



|V 2 f| 



(9) 



where || ■ || is the L2 norm on T 2 . It is apparent that 
we are restricting ourselves to sufficiently smooth forcing 
functions. 

The time and space averaged energy dissipation rate is 



,(|Vu| 2 ) = ,(c 2 ), 



(10) 



the second expression resulting from integrations by parts 
utilizing incomprcssibility. The bulk averaged enstrophy 
dissipation rate is 



X ^(|Va;| 2 )=^|V 2 u| 2 >- 



(11) 



We think of ft = e/kfU 3 and 7 = x/fcj£/ 3 as functions 
of Re and the functional form or "shape" of the forcing, 
but not explicitly on its amplitude 



F 



(12) 
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except indirectly through its influence on U . 

We are considering the Reynolds number to be the 
"control parameter" even though it is defined in terms 
of the emergent quantity U . Strictly speaking the flow is 
determined by the structure and amplitude of the body 
force (and the initial data) so the Grashof number such 
as Gr — F/kji/ 2 should naturally be used as the rele- 
vant dimensionlcss control parameter indicating the in- 
tensity of the driving and the resulting flow. Indeed, 
while we can always realize any given value of Gr, it is 
not at all evident that every particular value of Re can 
be achieved. Nevertheless, in order to express the results 
in terms of quantities familiar to the theory of homoge- 
neous isotropic turbulence we will (without loss of rigor) 
express the bounds in terms of Re. 

Poincare's inequality applied to HlOj) and I jlljl immedi- 
ately yield the lower estimates 



471 " tt2 



so that 



(3 > 4 7 r 2 a 2 i?e- 1 



and 



and 



X > v- 



16tt 4 



L 4 



■IT 



7 > 167r 4 a 4 i?e -1 



(13) 



(14) 



where a = (fc/L) -1 < (2tt)~ 1 is the ratio of the forcing 
to domain length scales. If (3 and 7 scale both as ~ 
Re^ 1 then we say that the flow exhibits laminar behavior 
because the energy is then necessarily concentrated at 
relatively long length scales determined by the prefactor, 
rather than over a broad range of scales that increases as 
Re — > 00. 

On the other hand if (3 ~ Re , the scaling expected 
in three dimensional turbulence, the flow exhibits finite 
(residual) dissipation in the limit of zero viscosity in- 
dicating the presence of an active and effective ene rgy 
cascade to small scales. It was recently shown [52, 
that (3 < cRe for the vanishing viscosity limit of three 
dimensional versions of the systems under consideration 
here and in section III, where the coefficient c is uni- 
form in is, L, and F. There is, however, no known a 
priori enstrophy dissipation rate bound for the three di- 
mensional turbulence; this is related to the outstanding 
open question of the regularity of solution for the three 
dimensional Navier-Stokes equations. As the results of 
this paper suggest quantitatively, the dissipation rates 
of two dimensional turbulence falls somewhere between 
laminar scalings and the rates for three dimensional tur- 
bulence. 

To prove the two dimensional bounds we first take the 
inner product of the vorticity equation © with to and 
average to obtain the enstrophy production-dissipation 
balance 1 

X = M) (15) 

where the time derivative term drops out when we take 
the long time average. Integrating by parts to move the 
k-Vx from to onto <f> and the Cauchy-Schwarz inequality, 
we easily obtain 

(16) 



X < kjUF. 



For the second step, consider a smooth incompressible 
vector field v(x, y). Take the inner product of v with the 
Navier-Stokes equation, integrate by parts and average 
to obtain 



1 

1? 



/ v ■ id 2 x = - (u • (Vv) ■ u + uu ■ V 2 v) . (17) 
Ji 2 , 



Using the Cauchy-Schwarz and Holder's inequality (as in 
p^ l we deduce 



F x / v • f d 2 x < C/ 2 ||Vv| 

L ll f ll Jn 



V 



|||V 2 v|| (18) 



norm on T 2 L . In order for the 



where || • is the L 
inequality to be non-trivial we need to restrict v so that 
J T 2 v ■ f d 2 x > 0. This is easy to arrange. For exam- 
ple the choice v = i/F will satisfy this condition if f is 
sufficiently smooth that the right hand side of (|18fl is fi- 
nite. If it is not so smooth, then for instance we can take 
v ~ A" o f where K(x, y, x' , y') is a (positive) smoothing 
kernel. In any case we can choose v appropriately and 
use l|18|l to eliminate F in Q16|l so that: 



X < U 3 k 3 f 



Ci + 



9i 

Re 



7< Ci 



Re 



(19) 

where the dimensionless coefficients Ci and C2 are inde- 
pendent of kf and L, depending only on the functional 
"shape" of v (and thus also on the shape of f ) but not 
on its amplitude F or the viscosity v. Explicitly they are 



C x = 



iVivll 



; v -f/F) 



and 



Co = 



|vH> 1/2 (20 ) 



where V; is the gradient with respect to the non- 
dimensional coordinate fe/X. An upper bound for the 
enstroph y d issipation rate like that in Ijl9(l was first de- 
rived in [26j. Note that for strictly band-limited forces, 
i.e., if the Fourier transform of the force is supported on 



wavenumbers with |k| G (fc. 



n ) k max ) 



with < k r 



< 



k m ax < °Oj then the coefficients G\ and Ci are bounded 
by pure numbers depending only on k max /k m i n . 

For the final step of the proof we use integrations by 
parts and the Cauchy-Schwarz inequality to see that 



,2\2 



= (u- V x (few)) < (|u| 2 )(|Vw| 



Combining (|2T|) with (fTi^l we deduce 

(u; 2 ) 2 < 



2,2 ^ k f U 5 („ , C 2 



Re 



(21) 



(22) 



and in terms of the energy dissipation rate this is the 
announced result 



,< Re -v*( Cl + g) 



1/2 



(23) 
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III. TIME-DEPENDENT FORCES 



Then using this to eliminate F from (|28[) we have 



Now consider the Navier-Stokes equation j£| where the 
time dependent body force f (x, y, t) is smooth and incom- 
pressible with characteristic length scale ~ kj 1 given by 



1.4 = 



iv 2 fr 



and time scale fij defined by 



We define 



(\dtf\ 2 
(Ifl 2 ) 



k t U' 



(24) 



(25) 



(26) 



the ratio of the "eddy turnover" time (kfU) 1 to the 
forcing time scale ■ In this time-dependent setting 
the amplitude F of the force is 



F= (Ifl 2 ) 1 / 2 . 



(27) 



As before, the space and time average of u> times 
the vorticity equation JfjJ yields the enstrophy balance 
equation (|15fl . and integration by parts and the Cauchy- 
Schwarz inequality implies 

X < kjUF. (28) 

For the second step here we introduce a smooth incom- 
pressible vector field v(x, y, t) and take space and time 
average of the inner product of with the Navier-Stokes 
equation to obtain 

(v f) = - (u • d t v + u • (Vv) • u + vu • V 2 v) . (29) 

Cauchy-Schwarz and Holder's inequalities then imply 

•f) 



F J^Iy75 < U(\d t v{ 2 y/ 2 + f/ 2 sup||Vv|| 
+ isU(\V 2 v\ 2 ) 1/2 - 



(30) 



Now we need to be able to choose v satisfying (v ■ f) > 
such that all the coefficients on the right hand side are 
all finite. Our ability to do this depends on details of 
f(x,y,t). 

For example if f is sufficiently smooth in space and 
appropriately uniformly bounded in time then we can 
choose v ~ f . We could also choose v as an appropriately 
filtered version of f to cover more general cases. For the 
purposes of this study and to display the results in the 
clearest (if not the sharpest or most general form) we will 
simply presume that f is sufficiently regular that we can 
take v = f . In that case (|30|l becomes 



F < tt f U + U 



2 SU Pt 



iVfl 



vk 2 U. 



(31) 



c kp 3 ( r + C 3 + i- ) 



i 

Re 



where the coefficient C3 is 



C, = 



F 



(32) 



(33) 



with V; denoting the gradient with respect to the non- 
dimensional coordinate fc/x. The dimensionlcss number 
C3 is independent of the scales of F, kf, L, etc., depend- 
ing only on the "shape" of f . For example if f is quasi- 
periodic with N frequencies and involves only wavenum- 
bers k with < k m i n < |k| < k max < 00, then C3 is 
bounded by V7V times a function of k max /k m i n . 
The final step again uses the inequality 

e 2 = v 2 {uj 2 ) 2 = v 2 {u ■ V x (kco)) 2 < vU 2 X (34) 

and it then follows immediately from l|32() that 

1/2 



0<Re~^ (r + C 3 + - 



(35) 



Note in this case r depends on U and features of the 
forcing through kf and 0/, but not on v. 



IV. MONOCHROMATIC AND CONSTANT 
FLUX FORCES 

An even sharper scaling bound on the energy and en- 
strophy dissipation rates can be derived when the driving 
is monochromatic in space, whether it is steady or time 
dependent 0, |53 ■ Suppose the body force involves only 
a single length scale, i.e., 



-V 2 f = k)i . 



(36) 



This docs not preclude complex timc-dcpcndcncc for 
f(x,y,t), just that it involves only spatial modes with 
wavenumbers k with |k| = kf. Then the enstrophy 
production-dissipation balance (|15fl implies 

X = {^) = (u • (- V 2 f )) = k) (u • f) = k 2 e. (37) 

Combining this with Q34[). we observe that 



e 2 < vU 2 x = vk)U 2 e (38) 



so that 



e < vk 2 f\J 2 



and 



X < vk)V 2 



(39) 



implying that both j3 and 7 are bounded by Re . Note 
that this kind of monochromatic forcing is a special case 
that has been shown in the literature for some cases to 
lead to a laminar flow that never looses stability [3^ ■ 
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An other type of forcing that results in this scaling is 

Vu 



f{x,y,t) 



L~ 2 Vu 



(40) 



where V is the projector onto spatial modes of wavenum- 
ber k with |k| <= [k m i n , k max ], and the coefficient e is now 
the control parameter. This type of forcing is often ap- 
plied in numerical simulations of homogeneous isotropic 
turbulence. With this forcing in the Navier-Stokes equa- 
tions constitutes an autonomous dynamical system with 
kinetic energy injected at a constant rate e at wavenum- 
bers with |k| <S [k m i n ,k max \. The rms speed U (i.e., the 
Reynolds number) and the enstrophy dissipation rate \ 
are then emergent quantities, determined by the imposed 
energy flux e. The mean power balance for solutions is 
still 

u(\Vu\ 2 ) = u(lu 2 ) = e, (41) 
and the enstrophy production-dissipation balance reads 

iwuir 



\Vu\\ 



(42) 



Because the the forcing only involves wavenumbers 
in [k m i n , k m ax] with positive energy injection at each 
wavenumber, at each instant of time 

kLnM 2 < l|VPu!| 2 < k 2 rnax \\Vu\\ 2 . (43) 

Then implies that 



h 2 i= < -v- < h 2 e 

Using this with (j34|) we see that 

e 2 < vU 2 X < vU 2 k 2 max e, 
and we conclude 



e < vk 2 n JJ 2 



and 



X < vk max U . 



(44) 
(45) 
(46) 



Hence also in this case both (3 and 7 are bounded ~ Re -1 . 

Note that in both these derivations a condition like 
(|44|l or the stronger condition i|37|) was used. It is an 
open question whether such a condition holds for more 
general and more "realistic" forcing functions. The re- 
sults <|32[1 and l|35|) give restrictions on the energy and 
enstrophy dissipation rate for a broader class of driving, 
but it is natural to wonder how broad of a class of forc- 
ing functions would actually result in the i?e _1 scalings 
in the vanishing viscosity limit. 



DISCUSSION 



U, L, kf and Qf. To be precise, e vanishes at least as 
fast as i?e -1 / 2 in this limit. This confirms that there is 
no forward energy cascade in the steady state in the in- 
viscid limit. On the other hand the residual enstrophy 
dissipation allowed by l|32|) in this limit docs not rule out 
a forward enstrophy cascade. This combination, e — > 
with \ = 0(1) in the inviscid limit, is consistent with the 
dual-cascade picture of two-dimensional turbulence de- 
veloped by Kraichnan 0, Leith an d Batchelor [Tlj . 
Note however that the absence of forward energy cas- 
cade is not necessarily true for any finite value of Re. 
The (3 ~ Re^ 1 / 2 scaling allowed by the bound is less 
severe than what a laminar flow (or a flow with only in- 
verse cascade of energy) would predict, and as a result 
(|23[l does not exclude the presence of some direct sub- 
dominant cascade of energy when the Reynolds number 
is finite as suggested by \u\ . 

On the other hand the direct cascade of enstrophy 
is necessarily absent for some forcing functions (see 
[l6l I27L l28j ) . When the forcing acts at a single scale or 
constant power is injected in a finite band of wavenum- 
bers, both e and \ vanish ~ Re^ 1 . This suggests an 
essentially laminar behavior for these flows: if the en- 
ergy spectrum follows a power law E(k) ~ k~ a for large 
wavenumbers then the exponent must be a < — 5 for \ 
to vanish in the vanishing viscosity limit. These results 
have been interpreted as absence of enstrophy cascade 
in finite domains. However, both of these results rely 
on the condition l|44|) which is not guaranteed for a gen- 
eral forcing functions. Whether 144fl might hold for more 
general forcing functions is an open question; the results 
(|32fl and l|35|) give the restrictions on the energy and en- 
strophy dissipation rate for a general forcing. Note that 
the (3 ~ Rer 1 ! 2 does not impose significant restriction 
on the energy spectrum given the bound on \. These 
considerations suggest that it is possible therefore that in 
two dimensional turbulence the steady state energy spec- 
trum depends on the type of forcing used, even within the 
class of relatively narrow-band driving. High resolution 
numerical simulations with forcing that does not neces- 
sarily satisfy (|44f> would be useful at this point to resolve 
this issue. 

We conclude by noting that an interesting question 
that follows from these results is that of the -Re-scaling 
of the energy dissipation in systems that almost have 
two dimensional behavior like strongly rotating, strongly 
stratified or conducting fluids in the presence of a strong 
magnetic field. For example, is there a critical value of 
the rotation such that the scaling of the energy dissi- 
pation rate with the Reynolds number transitions from 
e ~ Re to £ ~ i?e -1 / 2 ? These questions remain for 
future studies. 



These quantitative bounds show that for two dimen- 
sional turbulence sustained by forces as described in the 
previous sections, there is no residual dissipation in the 
vanishing viscosity limit defined by Re — > 00 at fixed 
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